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ABSTRACT 

We analyse cosmological hydrodynamic simulations that include theoretically and ob- 
servationally motivated prescriptions for galactic outflows. If these simulated winds 
accurately represent winds in the real Universe, then material previously ejected in 
winds provides the dominant source of gas infall for new star formation at redshifts 
z < 1. This recycled wind accretion, or wind mode, provides a third physically distinct 
accretion channel in addition to the "hot" and "cold" modes emphasised in recent the- 
oretical studies. The recycling time of wind material (ircc) is shorter in higher-mass 
systems owing to the interaction between outflows and the increasingly higher gas 
densities in and around higher-mass halos. This differential recycling plays a central 
role in shaping the present-day galaxy stellar mass function (GSMF), because declin- 
ing tree leads to increasing wind mode galaxy growth in more massive halos. For the 
three feedback models explored, wind mode dominates above a threshold mass that 
primarily depends on wind velocity; the shape of the GSMF therefore can be directly 
traced back to the feedback prescription used. If we remove all particles that were 
ever ejected in a wind, then the predicted GSMFs are much steeper than observed. 
In this case galaxy masses are suppressed both by the ejection of gas from galaxies 
and by the hydrodynamic heating of their surroundings, which reduces subsequent 
infall. With wind recycling included, the simulation that incorporates our favoured 
momentum-driven wind scalings reproduces the observed GSMF for stellar masses 
10^ Mq ^ M ^ 5 X 10^° Mq. At higher masses, wind recycling leads to excessive 
galaxy masses and star formation rates relative to observations. In these massive sys- 
tems, some quenching mechanism must suppress not only the direct accretion from 
the primordial IGM but the re-accretion of gas ejected from star-forming galaxies. In 
short, as has long been anticipated, the form of the GSMF is governed by outflows; the 
unexpected twist here for our simulated winds is that it is not primarily the ejection 
of material but how the ejected material is re-accreted that governs the GSMF. 

Key words: galaxies: formation, assembly, low redshift; intergalactic medium; cos- 
mology: theory; methods: numerical 



1 INTRODUCTION 

A fundamental challenge in galaxy formation theory is to 
understand the mechanisms that shape the observed galaxy 
stellar mass function (GSMF). This requires understanding 
how galaxies form and grow within dark matter halos. Ob- 
servations and theory indicate that the relation of galactic 
stellar components to underlying dark matter halos is not 



simple. In the mass range corresponding to halos of individ- 
ual galaxies, the dark matter halo mass function, $(Mh), 
follo ws roughly dn/dMh ^ (e.g ISheth fc TormenI 

1999| : IJenkins et al]|200ll : ISpringel et al.ll2005l ). turning over 
only in the group mass range. On the other hand, the 
GSMF, $(Ms), is most often fit using a Schech ter function, 
dn/dM, ~ M,"exp(-Ms/M*) (jSchechteJIlOT^ '). where Ms 
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is stellar mass and 2 = M* is roughly the mass of the 
Milky Way. The low-mass slope a is much shallower than 
-2 (e.g.'Bal drv et a\\ \200§). Hence the efficiency of star for- 
mation decreases at masses both lower and higher than M* . 
A theory of galaxy formation must explain the following fea- 
tures of the GSMF: a preferred mass for star formation ef- 
ficiency, a shallow low-mass end slope, and a rapid drop at 
large masses. 

We categorise physical mechanisms that imprint fea- 
tures on the GSMF into two classes: those having to do 
with how galaxies obtain baryons ("accretion"), and those 
having to do with how galaxies lose baryons ("feedback"). 
For accretion, the cooling of the hot virialised atmospheres, 
i.e. "hot mode" accretion, was historically believed to be 
the main channel of gas supply to galaxies. Models of this 
process give a characteristic mass above which gas in a 
gravitationally contracti ng halo cannot cool in a Hubble 
time (tu) to form s tars (|Rees fc Ostrikeij[l977l : [siii3ll977l : 
I White fc Re"3ll978l '). thereby providing a turnover in the 
theoretical GSMF, though in the absence of feedback the 
implied m ass scale is h igher and the cutoff shallowe r than 
observed jWhite fc Fren k 1991; iThoul fc Weinberd 1 19951 : 
iBenson et al.ll2003l ). 

More recently, "cold mode" accretion was argued to 
be responsible for most of the gas supply into g alaxies 
across cos mic time (IKatz et al. 2003; Keres et a l. 2005, here- 
after K05; iDekel fc Birnboim 2006). In cold mode, gas is 
supplied via cold filaments that rapidly stream gas from 
the intergalactic medium (IGM) into star- fo rming regions , 
without ever encounter ing a virial shock (|Binnevl 1 19771 : 
iBirnboim fc Dekell [200I) . This is most effective in lower- 
mass halos, while in halos > 10^^ M0 shocked gas domi- 
nates the halo gas, although dense filaments at high redshifts 
{z > 2) still allow cold mode accretion to proceed even above 
this threshold mass (K05; Dekel fc Birnboim 2006; Keres 
2007; Ocvirck et al. 2008; Keres et al. 2009a, hereafter K09a; 
Brooks et al. 2009; Dekel et al. 2009). At lower redshifts, the 
supply of cold gas declines in massive halos, and the cooling 
of the hot atmosphere becomes relatively more important 
(K09a). This late hot mode is relatively inefficient and can 
be naturally prevented in the large, constant density cores 
of cluster halos. 

Yet accretion alone cannot be the entire story, as it 
is well known that without strong feedback, models fail to 
produce basic ob servations such as the luminosity function 
and the GSMF (iKauffmann et al.l Il993l: ICole et al.l Il994 : 



iMuraU et all 12002 : iBower et al.ll2006l l. iKeres et all (|2009b . 
hereafter K09b) explored a cosmological SPH simulation 
without any form of explicit feedback and derived a "cor- 
rection factor" as a function of stellar mass to correct the 
simulated GSMF to match the observed one. The fact that 
their correction factor is always greater than one exempli- 
fies the "over-cooling" problem in galaxy formation, where 
the efficiency of forming stars, based on simple theoret- 
ical expectations, is greater than that act ually observed 
fe.g. lWhite fc Frenkll99ll : lBalogh et al.ll200ll 'l. Distinct feed- 
back processes may operate at different mass scales, given 
that star formation needs to be suppressed by ~ 10 at 
Ms < 10^° M0, falling to only ~ 2 at lO" M©, and then 
increasing again to ~ 5 by 10^^ Mq (see /corr plotted in 
K09b, Figure 1). 

K09b differentiated "preventive" feedback, where a pro- 



cess prevents gas from entering a galaxy's ISM in the first 
place, and "ejective" feedback, where gas is accreted and 
then thrown back out. Preventive feedback is most applica- 
ble for hot mode, since halos with dilute hot atmospheres 
can be more effectively kept hot by feedback processes such 
as jets from act ive galactic n uclei (AGN) and magnetic 
conduction (e.g. iBinne vl |2004 K05). An example of pre- 
ventive feedback affecti ng hot mode is "r adio" -mode feed- 
back (|Best et a l. 2005; Croton et al]|2006l '). However, K09b 
showed that preventing hot mode accretion alone is insuffi- 
cient to correct the GSMF, since in large systems most of 
the stars are still formed via cold mode when the galaxy was 
smaller (K09a), meaning that cold mode must be suppressed. 
Preventing cold mode accretion is much more challenging, 
since this gas is streaming in via dense filaments that are 
not easily disrupted. Hence, even for galaxies that are today 
dominated by hot mode accretion, ejective feedback appears 
necessary to understand their present-day stellar masses. 

In short, both the high and low ends of the observed 
GSMF appear to require an ejective feedback mechanism 
to curtail cold mode-driven star formation. A leading 
candidate for suc h a mechanism is galactic winds (e.g. 
iDekel fc SilSligsd ). Winds ejecting galactic ISM have the 
added benefit of chemically enriching the IGM, where 
metal absorption lines are observed to be common (e.g . 



Songailal '2001: Sc have et al] 120031 : 
Simcoe et al.ll2004 : .Adelberger et al. ,.,2005: ') . Star formation 
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suppression, IGM enrichment, plus the ubiquity of escaping 
outfiows o bserved in star forming galaxy spectra (e.g., 
Pettini et al. 2001; Shaplcv~eralT [2003: Martin 200M; 



Weiner et al.l |2009j) have motivated the inclusion of star 



formation bas ed feedback prescrip ti ons in cosmological 



simu lati ons ijSpringel fc Hernguistl 2003bl: Kobavash| 



2004 ICen fc Ostriker 2006; 'Oppenheimer fc Dav| 
20061: ISommer-Larson fc Fynbo 2008; Wicrsma et al] 
2009bl ). lOppenheimer fc Pavel (^06l ) explored a 
range of feedback prescriptions, finding that the scal- 
ings predicted by a momentum- driven wind model 
(|Murrav. Quatert, fc Thompsonl|2005l ) were most promising 
for matching the observations of high-2 IGM enrichment. 
In this scenario, wind velocity (iiwind) scales with galaxy 
velocity dispersion (a), leading to the gradual enrichment 
of the IGM matching the observed frequenc y of 2; ~ 6 C iv 
lines (|Oppenheimer. Dave, fc Finl ator 200 9|), the evolutio n 
of Civ from z = 5 -> 2 (|Oppenheimer fc Davi 120061 ) . 
plus the kinema tics of damped Lya absorbers (DLAs, 
iHong et al.ll2009l . in preparation). The mass loading factor 
{rj, where Mwind = vMsf) scales as a"^, meaning that 
lower mass galaxies are the most prolific enrichers by mass. 
This relation is central to reproducing the z •--^ 6 galaxy 
lumin osity function (|Dave. Finlator. fc Oppenheimeii 
l2006f), explaining the galaxy mass-metallicity relation 
(Finlator & Dave 2008), and providing a solution to the 
missing halo baryon problem (|Davell2009l ). The prescription 
even appears important for systems in the local Universe, 
as these wind simulations can accou nt for the 2 = missing 
metals (iDave fc Oppcnheimcr"2007^. the observed Ovi for- 
est at 2 < 0.5 (Oppenheimer fc Dave 200 9), and intragrou p 
enrichment and entropy levels at 2 = (jPave et alJbOOsI ). 
While our simulations are only sensitive to the predicted 
wind scalings, and do not try to model the detailed physics 
of wind propagation, there is independent evidence that 
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other fundamental properties of galaxies a re best under- 
stood by invoking momentum-driven winds (|Hopk ins et alj 
l20ld ). 

One consequence of having highly mass-loaded, 
moderate-velocity outflows emanating from lower mass 
galaxies is that wind material remains in the relative prox- 
imity of their parent galaxies. lOppenheimer fc Pavel (|2008l . 
hereafter OD08) found that even when outflows launch in 
excess of a halo's escape velocity, wind material regularly 
"recycles" back into a galaxy on a median timescale of 1-2 
Gyrs. This occurs because winds are more slowed hydrody- 
namically by interactions with halo gas and gas outside ha- 
los rather than gravitationally. Winds typically reach ~ 100 
physical kpc at all epochs, which means that outflows from 
^ AI* galaxies do not escape from their halo after 2 ~ 1. 
OD08 introduced the term halo fountain to describe such 
winds that remain trapped in their parent halos (see also 
iBertone et aDl2007l ). K09b coined mtergalacUc fountain to 
describe the transfer of feedback gas from low-mass galaxies 
at early times to higher mass galaxies at later times via a 
journey through the IGM. 

A major implication of our wind models is that the 
amount of material ejected from galaxies is, globally over 
cosmic time, far larger than the amount of material forming 
into stars. Hence these outflow fountains recycling back into 
galaxies can in principle represent a significant reservoir for 
stellar mass formation over a Hubble timescale. This paper 
is devoted to examining such accretion, which we call "recy- 
cled wind accretion" (or just "wind mode"), in cosmological 
hydrodynamic simulations. In addition to hot and cold ac- 
cretion, wind accretion effectively provides a third avenue for 
how galaxies get their gas. In this work we examine issues 
such as: How much of the cosmic stellar mass growth can 
this accretion mode create? Can this mode explain the pre- 
ferred mass of galactic star formation at M* (i.e. the knee in 
the Schechter function) , as well as the faint-end slope of the 
GSMF? What happens to the GSMF if winds never recycle? 
Does wind mode need to be "quenched" at any particular 
mass, and if so what physics might produce this quenching? 

In §2 we introduce our latest cosmological hydrody- 
namic simulations with several different prescriptions of su- 
perwind feedback. We begin to answer some the above ques- 
tions in §3 when we consider stellar mass growth divided into 
these component modes. We emphasise the concept of dif- 
ferential wind recycling as a function of galaxy mass and 
link it to the mass dependence of wind mode accretion. The 
effects of wind recycling on the GSMF are analysed in i)3.4l 
Our purpose is not to precisely fit the z — GSMF, but to 
introduce the concept of wind mode as a source of galaxy 
growth. We do argue that plausible wind models can come 
close to matching the observed GSMF, but only for galax- 
ies below M* . In §4 we discuss the balance of ejective and 
preventive feedback with recycled wind accretion. Ignoring 
wind mode, we examine how feedback suppresses star for- 
mation in finding that no feedback implementation we 
explore can adequately explain the observed GSMF with- 
out wind mode. We then emphasise the effect of wind mode 
on the slope of the GSMF in ^4.21 ^4.31 examines whether 
wind mode is partly a numerical effect, especially above M* 
where star formation-driven outflows are completely ineffec- 
tive at quenching star formation. Finally, §5 summaries our 
findings. 



2 SIMULATIONS 

We employ the Gadget-2 N-body + Smo othed Particle Hy- 
drodynamic (SPH) code (jSpringell 120051 ) to evolve a series 
of cosmological simulations to z = 0. We adopt a ACDM 
cosmology using cosm ological parameters b ased on the 5- 
year WMAP results (|Hinshaw et al.l \2QQ^ : = 0.28, 

r^A = 0.72, h = i^o/(100 kms"^Mpc"^) = 0.7, a primordial 
power spectrum index n — 0.96, an amplitude of the mass 
fiuctuations scaled to erg = 0.82, and = 0.046. We refer to 
these simulations as the r-series, where our general naming 
convention for a simulation is r[boxstze]n[particles/side][wind 
model]. Our primary simulations use a boxsize of 48 h~^Mpc 
with 384 dark and SPH particles/side resulting in a dark 
matter particle mass of 1.8 x 10* Mq and an SPH parti- 
cle mass of 3.6 x 10^ M©. Gas particles on average spawn 
two star particles resulting in star particle masses averaging 
1.8 X 10^ Mq. We define the wind model suffixes at the end 
of this section. 

An overview of the Gadget-2 code can be found in 
§2 of K09a. Additions to the public version of the code 
include cooling p rocesses usi n g the primordial abundances 
as described in iKatz et al.l (Il996|) and metal -line cool- 
ing as described in Oppenheimer fc Pavel (|2006l ). Star for- 
mation is modelled using a subgrid recipe introduced in 
ISpringel fc Hernguistl l|2003al ) where a gas particle above a 
critical density is modelled as a fraction of cold clouds em- 
bedde d in a warm ionised medium as in iMcKee fc Ostriker 
I1977I). The star formation follows a Schmidt law (Schmidt 
[1959I) with the star fo rmation timesca le scaled to match the 
z = Kennicutt law (jKennicuttll 19981 ) . The density thresh- 
old for star formation is nn = 0.13 cm~'^. 

We use a IChabried l|2003l ) initial mass function (IMF) 
throug hout, whi c h has a turnover at masses < 1 Mq relative 
to the ISalpeteJ (|l955h IMF. The fraction of mass in the 
IMF going into massive stars, defined here as ^ 10 Mq and 
assumed to result in Type II supernovae (SNe), is /sn ~ 
0.18 in the Chabrier I MF. As our /sn is larger than for the 
Salpeter IMF used in IS pringel fc Hernguistl (|2003al . /sn = 
0.10), the amount of SN feedback energy coupled to the ISM 
is greater per unit star formation as is the energy input into 
galactic outflows in our standard constant-wind model (see 
below). 

Our simulations directly account for metal enrichment 
from sources including Type II SNe, Type la SNe, and AGB 
stars. Gas particles eligible for star formation undergo self- 
enrichment from Type II supernovae (SNe) using the instan- 
taneous recycling approximation, where mass, energy, and 
metallicity are assumed to immediately return to the ISM. 
Type II SN metal enrichmen t uses the metallicit y -depe ndent 
yields calculated from the IChieffl fc Limongil (|200i) SNe 
models. Our prescriptions for feedback, described below, also 
assume the instantaneous input of energy from O and B 
stars (i.e. sta rs that are part of /sn). We inp ut the Type la 
SNe rates of IScannapieco fc BildstenI (|2005h . where an in- 
stantaneous component is tied to the SFR, and a delayed 
component is proportional to the stellar mass as described 
in OD08. Each Type la SN adds 1 0^^ ergs of thermal energ; 
and the calculated metal yields of iThielemann et al.l (|l98' 
to surrounding gas, while the mass returned is compara- 
tively negligible. AGB stellar enrichment occurs on delayed 
timescales from 30 Myrs to 14 Gyrs, using a star particle's 
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age and metallicity to calculate mass loss rates and metallic- 
ity yields as described in OD08. AGB mass and metal loss 
is returned to the three nearest surrounding gas particles. 
OD08 showed that the largest effect of AGB stars is to re- 
plenish the ISM, because a star particle can lose nearly 40% 
of its mass over tn, about double the /sn ~ 18% that is 
recycled instantaneously via Type II SNe. 

Feedback is directly tied to the SFR, using the relation 
Mwind = ?7MsF, where rj is the "mass loading factor." The 
probabilities for an SPH particle with a non-zero SFR to 
turn into a star particle and to join a wind are calculated in 
each timestep; the ratio of the two probabilities is therefore 
rj. Actual conversion into a star or wind particle is deter- 
mined stochastically using these probabilities and a random 
number generator. We explore simulations with no winds 
and three different forms of implemented feedback. 

No Winds (nw): No feedback energy is imparted ki- 
netically to SPH particles. However, as with all simulations, 
ene rgy is imparted thermally to ISM SPH particles using 
the lSpringel fc Hernauis"^ l|2003al ) subgrid two-phase recipe 
where all SN energy is instantaneously coupled to the ISM. 
The r48n384nw simulation can be compared to the Gadget- 
2 no wind simulation explored in K09a and K09b, although 
there are some notable differences. First, the cosmology is 
different; they use = 0.26, = 0.74, h = 0.71, n = 1.0, 
as = 0.75, and ilb = 0.044. Second, the r-series uses a 
Chabrier IMF with delayed mass loss from AGB stars and 
feedback from Type la SNe, while the K09a simulation uses 
a Salpeter IMF and only Type II SNe feedback. The recy- 
cling of gas from AGB stars causes more late-time star for- 
mation in the r-series. A third difference is the addition of 
metal-line cooling for non-star forming SPH particles. This 
actually has little effect in the no wind case, because a vast 
majority of metals 96%) remain within galaxies in ei- 
ther the ISM or star particles. We will find throughout §3 
that the r48n384nw simulation reproduces the main mass 
growth properties (i.e. cold & hot mode) of the 50 /i~^Mpc, 
2 X 288^ particle simulation of K09a/b, although with some 
minor differences. 

Constant Winds (cw): Feedback energy is added ki- 
netically to gas particles at a constant rate relative to the 
SFR, with ?7 = 2 and a constant velocity Hwind = 680 kms~^. 
This is intended to mimic the constant-wind model of 
ISpringel &: Herngui st' ('2003a'), however with some notable 
differences. Here ri is defined relative to the SFR of the entire 
Chabrier IMF, and not j ust the long-lived stars in a Sa lpeter 
IMF as it is defined bv lSpringel fc Hernguist l|2003al '). Sec- 
ond, Uwind ~ 680 kms~^ is used rather than 484 kms~^, be- 
cause there is more SN energy per mass formed in a Chabrier 
IMF. The kinetic energy imparted to wind particles per unit 
star formation is 9.25 x 10** erg/ Mq, which is 95% of the 
SN energy in this IMF assuming all stars ^ 10 Mq add 10^^ 
ergs/SN. Wind particles are hydrodynamically decoupled 
until either 2.9 x 10^ yrs has passed or, more often, the SPH 
part icle has reached 10% o f the st ar formation critical den- 
sity. ISpringel fc Hernguisj (|2003lJ ) demonstrated that this 
decoupling achieved resolution convergence. It is supposed 
to simulate the formation of hot chimneys extending out 
of a disk galaxy providing a low resistance avenue for SN 
feedback to escape into the IGM. 

Slow Winds (sw): We simulate an alternative 
constant-wind model with the only difference being the 



outflow velocities are half as high as the cw winds (i.e. 
'^wind = 340 kms~^ and rj — 2). Therefore, only a quarter of 
the SN energy couples kinetically into these outflows. 

Momentum-driven Winds (vzw): The 
momentum-driven wind model us es the scalings of 
iMurrav. Quatert. fc ThompsonI (|2005l ) based on the galaxy 
velocity dispersion (a), a nd matches the observed trends of 

^wind 

bvlMartinl (l2005al} . Its implementation is explained 
in lOppenheimer fc Davt (|2006 ) and updated in OD08 with 
the addition of an in-r un group fi nder to calculate a using 
galaxy mass following IMo et al] (|l998l V One modification 
from OD08 is the in-run group finder now uses dark matter 
in addition to baryons to calculate masses as this produces 
better agreement across time with our post-run group 
finder (discussed below). The wind parameters vary across 
galaxies using the following relations: 



^wind 

V 



= 3av//L - 1, 

— — 

a 



(1) 

(2) 



where /l is the luminosity factor in units of the galactic 
Eddington luminosity (i.e. the critical luminosity necessary 
to expel gas from the galaxy potential), and ao is the nor- 
malisation of the mass loading factor. Here we include a 
metallicity dependence additionally for /l owing to more 
UV photons being output by lower metallicity stellar popu- 
lations 



/l = /l,0 X 10 



-0.0029*{log -1-0. 417694 



(3) 



randomly select /l,0 = [1.05, 2] for each S PH particle, fol- 
lowing observations bv lRupke etlo] l|2005l ). and add an ex- 
tra kick of size 2a (the local escape velocity) to simulate 
continuous pumpin g of winds until it is far from the galax y 
disk (as argued bv lMurrav. Quatert. fc Thompson! l|2005h ). 
The only free parameter is ao = 150 kms~^, which is the 
same efficiency factor of the mass loading used in OD08. 
Similar to the cw model, hydrodynamic forces are decou- 
pled for a time of 1.95 x 10^" kms~^/wwind yr or until 10% 
of the star formation critical density is reached; the latter is 
the case for the majority of wind particles. 

The in-run group finder used to derive wind properties 
produces robustly consistent masses used to calculate a com- 
pared to the post-run halo finder, although there are differ- 

1/3 

ences. The velocity dispersion of a halo, (Jh, scales as AI-^ , 
where Mh is derived by post-run. ah is different than a calcu- 
lated by the in-run group finder, which uses a higher thresh- 
old density to identify individual galaxies in groups and clus- 
ters rather than the dispersion of the entire group/cluster 
halo. Therefore, we find simulated winds have median val- 
ues following 7? (X cr-°-^^ between Mh = 10^°'^ - 10^^'^ Mq. 
The Wwind medians scale as aJJ *^, deviating from a linear re- 
lation owing primarily to lower metallicities from low-mass 
galaxies resulting in an extra boost (Equation [3|; rj does 
not depend on this equation. Because Uwind does not lin- 
early scale with a^, winds launched from galaxies residing 
in groups and clusters are often launched below the halo es- 
cape velocity (wcsc); we explore this further in tj3.2l At z — 1, 
a lO^'^ M0 halo identified by the post-run halo finder pro- 
duces winds with median Uwind = 490 kms^^ and ry — 1.7, 
equivalent to 40% of the total SN energy. 

Our four simulations are listed in Table [T] They were 
run at the University of Arizona's SGI cluster, ICE, the Na- 
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Table 1. Simulations" 



Name ''^wind 

(kms-i) r,'' B^ind/^^g: 



r48n384nw - - - 

r48n384cw 680 2 0.95 

r48n384sw 340 2 0.24 

r48n384vzw oca oc (j-i 0.72 

" All simulations have 384"^ SPH and dark matter particles in 
a 48 comoving /i~^Mpc box, corresponding to 1.8 X 10* and 
3.6 X lO'^ Mq per SPH and dark matter particle, respectively. 
They were all evolved to 2 = 0. The equivalent Plummer 
gravitational softening length is 2.5 comoving /i~^kpc. 
' is defined as Mwind ocr} X Mgp . 

Feedback energy divided by SN energy. Averaged over resolved 
galaxies for the r48n384vzw simulation. 

tional Center for Supercomputing Applications' Intel clus- 
ter, Abe, and the University of Massachusetts' Opteron 
cluster. Eagle. Each simulation took between 100,000 and 
150,000 CPU hours, although nearly 1 million CPU hours 
were used when including smaller box sizes and various test 
simulations. 

We use the program to identify bound groups 

of cold baryons and stars (see K05 and K09a for more de- 
tails). We have modified SKID with the requirement that 
bound gas particles must have non-zero SFRs; therefore 
we only identify gas in the galactic ISM and not from 
the extended halo. Our galaxy stellar mass limit is de- 
fined as ^ 64 star particles, Ms ~ -^qS-os because 
iFinlator et al.l (|2006l ) showed stellar properties of galaxies 
are well- converged above this limit. Halos are identified us- 
ing a Spherical Overdensity algorithm with the same param- 
eters as in K09a. 



3 RESULTS 

3.1 Modes of Star Formation 

We begin our discussion of stellar mass growth by dividing 
gas accretion and star formation into its component modes. 
Just as in K05 & K09a/b, we divide hot and cold mode 
based on the maximum temperature, Tmax, an SPH particle 
achieves on its journey into a galaxy. The temperature split 
is the same - T = 2.5 X 10^ K, which K05 justified as a good 
empirical division of the two modes. We define a third mode, 
wind mode, as any particle previously ejected in a wind that 
is re-accreted by a galaxy and forms stars (i.e. recycled). 
There is no temperature constraint on wind mode. 

We have modified our version of Gadget-2 to track 
these three modes during a simulation using the following 
rules. An SPH particle's Tmax is tracked every timestep if it 
has not yet been launched in a wind. If the particle's current 
temperature is greater than Tmax, then Tmax is set to the 
current temperature. Tmax cannot be modified if a particle 
has a non-zero SFR, where t he temperature is an average o f 
the two-phase ISM model of ISpringel fc Hernguistl (|2003al ). 
Particle wind launches are also tracked, and Tmax is not 

^ http://www-hpcc.astro.washington.edu/tools/skid.html 



allowed to be adjusted after such an event. Any star forma- 
tion from a previously wind-launched particle is considered 
wind mode. For the purposes of this paper, we do not con- 
sider a particle's Tmax nor subsequent temperature history 
once it is launched in a wind. We also cannot properly track 
gas transferred via AGB feedback from star particles to sur- 
rounding gas particles, although this is a second-order effect 
that is not expected to favour any one mode; the transferred 
material will inherit the gas particle's mode. 

We consider the global star formation divided into these 
three modes for our four models in Figure[T] Our baseline for 
comparison is the no-wind (nw) simulation in the top panel. 
The black line is the global SFR density as a function of 
redshift, and the blue and red curves correspond to the hot 
and cold modes of star formation respectively; we use this 
colour scheme for these modes throughout this paper. The 
nw model is the focus of K09a and K09b, who use a lower res- 
olution Gadget-2 simulation without winds as mentioned 
in Sj2] Our nw simulation shows the same qualitative charac- 
teristics identified by K09a including: (i) cold mode accre- 
tion dominates at all epochs, and (ii) the hot mode accretion 
fraction increases at late times; we refer the reader to their 
§3 for an in-depth discussion. However, note that Figure [1] 
plots star formation and not accretion (as was done in Fig- 
ure 2 of K09a), although the two are tightly coupled (K05). 
Throughout this paper we consider only SFRs and not ac- 
cretion rates, as our focus is on stellar mass growth; we leave 
a direct investigation of accretion rates for future work. 

The next three panels of Figure [T] show the global SFR 
densities divided by mode for our three feedback models; the 
thin black line is the nw global SFR shown for comparison. 
In addition to the hot and cold modes indicated by the red 
and blue lines, the green lines in each of these panels indicate 
the contribution from the recycled wind mode. In every case, 
the wind mode dominates the global SFR over most of the 
history of the Universe. The wind mode contributes over 
50% of the star formation after z = 1.3, 2.7, and 1.7 for the 
cw, sw, and vzw models, respectively. The wind mode grows 
relative to the other modes at late times, making up 84%, 
90%, and 92% of the z = star formation, respectively. The 
resulting integrated stellar mass fractions in 2: = galaxies 
are 51%, 73%, and 71%, respectively. 

This leads us to our first main conclusion: wind material 
not only commonly recycles back into galaxies, it in fact 
dominates late-time stellar mass growth. Ejected material 
cannot be assumed to remain in the IGM forever. The wind 
mode and the total SFRs decline at late times, albeit at a 
slower rate than the cold and hot modes. SFRs are lower for 
every feedback simulation relative to the nw simulation, but 
the reduction is fairly minor for sw winds and vzw winds 
at later times. Outflows efficiently suppress both the cold 
and hot modes in all cases, but the wind mode re-supplies 
some of the ejected gas ba ck into galaxies. Compare d with 
the observed SFR density ijHopkins fc BeacomI 12006^) (grey 
dashed lines), these two models over-predict star formation 
after 2 > 1, which we will argue in H4.3l is primarily a result 
of excessive star formation in > M* galaxies. 

3.2 Differential Wind Recycling 

Wind material recycles back into galaxies in a manner that 
depends on galaxy mass (OD08). This trend, which we re- 
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Figure 1. The global SFR density ( Mq yr~i Mpc'^) in co- 
moving units as indicated by the black lines, and sub-divided by 
accretion mode: cold (blue), hot (red), and wind (green). The in- 
tegrated global no-wind SFR density is repeated with thin lines 
in the bottom three panels for comparison with each wind model. 
The grey dashed line indicates the Hopkins & Beacom. (2006) 
compilation calibrated to a Chabrier IMF. All galaxies (centrals 
and satellites) are included. 



fer to as differential recycling, is key to understanding how 
recycling shapes the galaxy stellar mass function in these 
simulations. Wind recycling is defined as an event where 
a previously launched wind particle is either launched in a 
subsequent wind or converted completely into a star parti- 
cle. The recycling time tree is defined as the time between the 
initial launch and re-launch/star particle conversion, which 
includes the time spent going out of the galaxy, the time 
spent coming back into the galaxy, and that spent within 
the new galaxy's ISM. OD08 showed that individual wind 
particles average 2-3 wind launches (i.e. 1-2 recyclings) in 



their vzw simulations with a median t^cc of < 1 Gyr, show- 
ing a weak dependence with redshift. 

OD08 also showed recycling in the vzw model was 
highly dependent on a galaxy's baryonic mass, with tree scal- 
mg as ~ h'C^ for vzw outflows. At face value this is sur- 
prising, since in this model winds are launched at speeds 
nearly proportional to Vcac of the galaxy. Hence if grav- 
ity were the dominant criterion for how far a wind travels 
and Wwind < Veac, ouc should expect almost no trend with 
mass. By comparing the trends with analytic expectations, 
OD08 argued that hydrodynamical slowing dominates the 
wind deceleration, not gravitational forces. The deceleration 
is greater for more massive galaxies because they reside in 
denser gas environments, leading to quicker recycling. OD08 
demonstrated that the gas outside halos is primarily respon- 
sible for hydrodynamic slowing in lower mass halos, while 
winds rarely escape the gas-rich halos of > M* galaxies at 
z <1 (see their Figures 13 and 14). 

Figure [5] plots tree as a function of the halo mass for 
central galaxies at the time they launch their winds. We 
focus on wind particles launched between z = 1 ^ 0.95, 
for several reasons: (i) most wind mode accretion occurs 
below z = 1, (ii) there is ample time afterward (~ 8 Gyrs) 
to follow the outflow material, and (iii) there is a larger 
dynamic range in A/h to explore dependencies than at higher 
2. The curves at other redshifts do not differ dramatically 
from 2 = 1 for the same mass galaxy (see OD08) . The black 
points correspond to individual wind particles; the vertical 
striping at higher halo masses results from the finite number 
of such halos. The un-recycled winds appear as the thin line 
of points at ~ 8 Gyrs. The lines show the times below which 
10, 25, 50, 75, and 90% of launched particles recycle. If a line 
exists at a given Mh, then at least its corresponding fraction 
of particles recycle by 2 = (e.g. if top dotted line exists, 
^ 90% of winds recycle for that Mh). Finally, the histograms 
along the bottom show the relative number of integrated 
wind launches from all halos within a Mh bin. The statistics 
here only include central galaxies, which contain 87% of the 
total 2 = 1 star formation in each of the three models. 

The purple curves show the analytical tree expected 
purely from gravity (i.e. freely falling trajectories) for 
each model. These are calculated by integrating the time- 
varying gravitational acceleration along a particle's trajec- 
tory using the fwind as the lau nch velocity and assuming 
iNavarro. Frenk. fc White! (|l997f) h alo p otentials with con- 
centrations from the Duffy et al.l (|2008l ) relation at 2 = 1 
and a launch radius of 0.0125 r2oo- tree is twice the time 
needed to reach the turnaround radius, and is plotted if a 
particle is calculated to recycle before 2 = 0; no evolution 
for the halo is assumed. The median ti-oc from the cw and 
sw models are shorter by x2 — 3 compared to the analyti- 
cal case owing to hydrodynamic forces when t^cc ^ 1 Gyr, 
confirming the argument of OD08 that hydrodynamic slow- 
ing plays a critical role in wind recycling. While winds in 
the vzw model are launched in excess of Vcsc at high-2, this 
is not necessarily the case for galaxies > A/* at 2 < 1. 
Winds launched with Uwind oc a, where a is calculated for 
an individual galaxy, cannot escape the group and cluster 
potentials in which such massive galaxies typically reside. 
We use the median «wind tracked in the simulation for a 
given 2 = 1 Mh, therefore there exists some slight scatter 
from the expected Uwind — Mh relation. This scatter becomes 
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Figure 2. Recycling times - the time between when a wind particle is launched and when it is converted to a star or launched for a 
second time - for all wind particles ejected between z = I — 0.95 in the three wind simulations, as a function of the host galaxy's at 
launch time. Black dots show individual wind particles, and coloured solid lines show the median tree; if this line exists for a given 
then at least 50% of wind particles recycle. Similarly, 25 and 75% cuts on tree are indicated by dashed lines, while 10 and 90% cuts are 
indicated by dotted lines; if the top dotted line exists, > 90% of winds recycle. Purple lines show the expected tree from gravitational 
considerations alone. Histograms show the relative number of wind particles per Mjj emanating from central galaxies. Wind events that 
escape until z = contribute to the thin line of points at ^ 8 Gyrs, the time between z = 1 — >■ 0. The bottom right panel shows the 
diversity in the differential (i.e., mass-dependent) recycling curves (medians). Dotted grey lines, applicable for all models, indicate the 
redshift when median wind particles, launched at (z) = 0.975, recycle. 



amplified when performing the calculation of gravitational 
tree, leading to the non-monotonic trend of the purple line in 
the upper right panel. Gravitational forces play a significant 
role in slowing vzw winds, yet hydrodynamic slowing still 
reduces free by x3 — 5. 

The shape of the frec-Afh relation for all the wind mod- 
els we examined can generally be described by a power law 
decay of free with mass followed by a gradual fiattening. 
The differences between models can be seen explicitly in 
the bottom-right panel of Figure [2] where the median t^cc 
curves are compared. The cw model exhibits a steeper de- 
cline relative to the vzw model, which is simple to under- 
stand. The cw model has no dependence of outflow velocity 
on mass, leading to a greater gravitational deceleration at 
higher masses as the wind velocity falls below the typical 
fesc of the halo. Another strong prediction of the cw model 
is, as expected, nearly no gas recycling occurs in halos with 
masses below Mh ~ 3 x 10^^ M0. This is in sharp contrast 
to the vzw model where even low-mass halos experience sig- 
nificant rapid gas recycling. The sw model shows similar, 
but appropriately shifted, behaviour as the cw model, where 
the tree — Mh relation has moved to lower mass correspond- 
ing to lower Vesc- Quantitatively, t^cc ~ M^^'^ in cw versus 
~ M^"'^ in vzw over the range where tree reduces from 5 to 
1 Gyr. Above Mh ~ 10^^ M©, the median tree in the vzw 
model flattens and never decreases much below 1 Gyr. This 
contrasts to OD08 where tree approaches 100 Myrs for the 
most massive galaxies. This difference owes to the fact that 



we impose no speed limit in the simulations presented here 
whereas OD08 limits the iSwind to be no more than 2i5sN, 
reducing the iiwind from massive galaxies. 

Finally, the histograms of wind particle launches show 
significant differences between the models. These histograms 
are tied to the integrated SFRs per halo mass bin in central 
galaxies using the conversion factor rj"^ . For the constant- 
wind models, this conversion factor is fiat by definition, al- 
lowing the histogram to be a proxy indicator of the inte- 
grated star formation in central galaxies. The up-tick in the 
histogram at higher masses corresponds to the steep reduc- 
tion in tree aud is the signature of wind recycling significantly 
increasing star formation in higher mass halos at z = 1. This 
up-tick occurs at lower Mh in the sw model relative to the 
cw model, which we will show is related to the Afh where 
the wind mode dominates in the next subsection. Also, the 
integrals of each histogram are comparable between panels, 
and are proxies for the global z = 1 SFRs in Figure [T] (not 
including the minor contribution from satellites) . The higher 
histogram of the sw model relative to the cw model results 
from the greater z — 1 wind mode SFRs. The vzw model 
cannot be so easily interpreted since ri declines with halo 
mass. 

To summarise, the wind velocity primarily sets the recy- 
cling time at a given halo mass, while the mass loading factor 
multiplied by the SFR determines the integrated amount of 
winds emanating from a given halo mass range (i.e. the his- 
tograms in Figure [2]). The wind launch histograms increase 
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sharply where tree becomes smaU relative to tu, which links 
differential recycling and the re- accretion of gas for star for- 
mation. 

3.3 Star Formation and Stellar Masses 

Having examined how winds recycle, we now consider this 
gas resupply as a source of star formation in Figure [S] where 
the SFR modes are decomposed as a function of Mh in our 
four models at three redshifts (z —3, 1, and 0). The data 
points represent individual galaxies, and the lines show the 
mean SFRs as a function of Mh. For the no- wind simulation 
in the top panels, we reproduce the main trends of K09a: 
(i) SFRs at a given redshift increase with halo mass, (ii) 
SFRs at a given Mh decline with cosmic time since at least 
z = 4, and (iii) cold mode star formation dominates at all 
redshifts and at almost every halo mass with the exception 
of intermediate masses at late times. 

The lower three rows in Figure[3]correspond to our out- 
flow simulations; all have wind mode SFRs that increase 
steeply with Mh (green lines), and which dominate above a 
certain Mh (except for the vzw model at z = 3). The key 
thing to note is that the wind mode begins to dominate at 
approximately the Mh where the median free equals the Hub- 
ble time (cf. Figure [2] and z — 1 SFRs). Above this mass, 
wind recycling provides a gas supply that quickly exceeds 
the feedback-suppressed cold and hot modes. 

Although outflows suppress the cold and hot modes in 
these simulations, it is the mass dependence of differential 
recycling that determines where and how much wind mode 
star formation occurs. For example, the cw SFR— Mh re- 
lations exhibit dramatic rises at Mh ~ 10^^ M© at all red- 
shifts, because the median free falls below tu above this mass. 
The sw SFRs increase steeply at Mh ~ 10^^ Mq for the same 
reason, resulting in the global SFR density becoming wind 
mode dominated at z — 2.7 compared to z = 1.3 for the cw 
(cf Figure [T]). In contrast, the vzw model does not show a 
preferred Mh where SFR dramatically increases, because the 
differential recycling curve is shallower and extends to lower 
masses. For the vzw model by 2 = 0, the wind mode dom- 
inates in halos above 2 x 10^^ M©, leading to more overall 
star formation in sub-M* halos relative to the cw model. 

The total SFRs of the wind models approach and even 
exceed the total no-wind SFRs (thin black lines) at higher 
masses a,t z — 1 and 0. In all our implementations, star 
formation-driven feedback appears unable to prevent star 
formation in massive halos. The sw SFRs are similar to the 
nw SFRs above Mh = 10^^ M©, at all redshifts. The sw 
winds, launched at 340 kms~^, produce halo fountains with 
tiec <S tu, which result in similar star formation rates as 
in the no-wind case. In some sense, this is s imilar to the 
canonical scenario for winds envisioned in e.g. iDekel fc Sili3 
l|l986ll , where star formation-driven feedback has a strong 
impact on low mass galaxies but little impact at high masses, 
albeit for different reasons. 

All the wind models show greater SFRs relative to the 
nw case in a range of halos by z = 0, with the sw and 
vzw SFRs exceeding the nw case for Mh > 10^^ Mq. In 
part, the enhanced late-time star formation reflects a de- 
lay, with wind ejection and recycling shifting star formation 
from early times to late times. In addition, wind recycling 
can shift star formation from low mass galaxies to high mass 



galaxies via intergalactic fountains. Either way, wind recy- 
cling provides a greater reservoir of gas for late-time star for- 
mation. The resulting wind mode-dominated SFRs in high 
mass halos far exceed observed SFRs in the local Universe. 

Figure |4] shows the cumulative contributions of cold, 
hot, and wind mode accretion to the z = stellar mass of 
central galaxies. For the nw simulation, our results are sim- 
ilar to K09a/b: cold mode dominates at every mass, with 
a median contribution near 100% at low masses, 75% at 
Mh ~ 3 X 10" Mq, and 80-90% above Mh ~ 5 x 10^^ Mq. 
While direct cold mode accretion dominates at low masses, 
the highest mass galaxies are built largely from mergers of 
low-mass galaxies that formed their stars from cold mode 
gas (K09a). To reproduce the observed (much lower) masses 
of these galaxies, their low-mass progenitors must efflciently 
suppress star formation at early times, likely using an ejec- 
tive feedback mechanism (K09b). 

In contrast to the nw case, in all three wind models it 
is the wind mode that dominates the growth of high mass 
galaxies. The onset of wind mode dominance is swift in the 
cw and sw models, and slower in the vzw model owing to its 
shallower differential recycling curve. The transition mass 
where wind mode accounts for over half of the stellar mass 
growth occurs at Mh = 10^^ ^ lO" ■^ and lO" * Mq for the 
cw, sw, and vzw simulations respectively. The correspond- 
ing transition masses for the z = SFR (right panels of 
Figure [3]) are lower because recycled wind star formation oc- 
curs predominantly at late times. We show the link between 
Mh and Ms for central galaxies in Figure |4] using vertical 
grey lines drawn at the median Mh for central galaxies of 
Ms = W^, W^°, 10", and 10^^ Mq. The grey lines move 
to higher halo masses for the 10^ and lO^'' Mq lines with 
the addition of feedback, because a galaxy with the same 
stellar mass will reside in a more massive halo owing to 
feedback suppression. The lines shift much less for 10^^ Mq 
central galaxies as wind recycling renders feedback ineffec- 
tive at limiting stellar mass growth for high-mass galaxies. 
In terms of stellar mass, wind mode accounts for over half 
of the stellar mass growth above Ms — 10^"'*, 10^°'^, and 
10'-* ** Mq for the cw, sw, and vzw simulations respectively. 

3.4 Galactic Stellar Mass Functions 

We now consider stellar masses for all galaxies. Figure [5] 
shows the simulated z = GSMFs including all accretion 
(thick lines) and excluding wind mode (thin lines). The lat- 
ter case is equivalent to assuming wind materials never recy- 
cle and remain outside the ISM of any galaxy. Because SPH 
particles are tagged when flrst launched, all subsequent star 
formation can be removed from our analysis if particles recy- 
cle more than on ce. For comparison, we show the observed 
iBell et all (|2003l ') GSME0- As in K09b, the nw simulation 
signiflcantly over-predicts the observed GSMF at all masses. 
The cw model changes from a nearly flat power law without 
wind mode (a = -2.11 between Ms = 10^ '^ - lO"'^ Mq) 
to having a bump above Ms = 10^^ Mq including all ac- 
cretion. The result is a poor fit everywhere. Adding wind 
mode in the sw model fiattens the GSMF the most between 

^ For alternative rec ent determinations o f the GSMF, see 
iLi fc Wliit j | |2009|'| and lBernardi et al.l | |2009|'| . 
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Figure 3. Instantaneous SFRs of central galaxies as a function of Afh for three redshifts (z = 3, 1, 0) in our four simulations. Black 
points show total SFRs of individual galaxies, and the black line shows the mean. Blue, red, and green lines show the mean SFR for 
cold, hot, and wind modes respectively. The thin black lines indicating the total nw SFR are repeated in the feedback model panels for 
comparison. 



Ms = lO''-^ - 10^°-^ Mq, from a = -2.18 to -1.11. The 
more gradual increase of wind mode with mass changes 
the vzw faint-end slope from a = —1.92 to —1.45 over 
Ms = 10^-^ - 10^°-^ Mq, and provides the best match to 
the observations at Ms ^ 5 x 10^" M©. Excluding wind 
mode yields acceptable agreement at the highest masses for 
the sw and vzw models, but it spoils the agreement at lower 
masses. We discuss the relative shapes of the various GSMFs 
further in g2] 

To better characterise the impact of winds and wind 
recycling on galaxy masses, we rank order the galaxies in 
each simulation by stellar mass, then compare the proper- 
ties of galaxies of the same rank in the different simulations. 



In a statistical sense, this procedure picks out the directly 
comparable galaxies between pairs of simulations. Figure [§] 
plots the suppression factor, /supp = Mnw/Mwind relative to 
the stellar mass in the no- wind simulation, Afnw. A value 
/supp = 1 (grey line in the figure) means that the rank- 
ordered masses between simulations with and without feed- 
back are identical, i.e. feedback has had no effect. The solid 
black curve shows the value of /supp required to obtain the 
mass of observed gala xies with the s ame comoving space 
density according to iBell et al.l l|2003f ) ; this is the quantity 
referred to as fcon in K09b (their figure 1). The no- wind sim- 
ulation produces galaxies that are too massive by a factor 
~ 2.5 at Mnw ~ 10^^ Mq, rising to a factor of 5 — 10 at the 
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els is considerably larger than the naively expected factor of 
1 + t; = 3 (dotted grey line) , for the reasons that will be dis- 
cussed in tj4.1| - the difference partly reflects the impact of 
mass recycling from evolved stars, but the dominant effect 
is that winds heat the surrounding gas and suppress accre- 
tion, providing preventive feedback in addition to ejective 
feedback. The suppression in these models is relatively in- 
dependent of mass, compared to the steadily dropping /supp 
in the vzw model owing to the rj oc scaling of the wind 
mass loading. 

Including recycled wind mode (as our simulations do 
by default) changes the picture radically, as shown by the 
thick curves in Figure |6l For the cw and sw models, wind 
recycling has no impact on galaxies with masses below a 
threshold Mnw For the slow winds of the sw model, /supp 
approaches one above Mnw ~ lO" Mq. The cw model shows 
a similar behaviour at an ~ 8 x higher mass, but these faster 
winds suppress the masses of even the largest galaxies by a 
factor of two. The transition between the low and high mass 
regimes occurs at approximately the scale where ticc ~ tn in 
Fi gure[^ /supp in the vzw model also approaches one above 
Mnw ~ 2 X 10^^ Mq, but it continues to rise all the way to 
small masses unlike the constant wind models, which flatten 
towards the smallest masses. 

The full /supp curves show a much stronger mass de- 
pendence and much larger differences from model to model 
than the curves without wind mode. We conclude that it is 
the mass dependence of the recycled wind mode, what we 
have called differential recycling, that plays the strongest 
role in determining the shape of the GSMF in our simula- 
tions. In all the models, wind recycling leads to excessive 
stellar masses and star formation rates in massive galaxies 
compared to those observed. However, simply suppressing 
all wind recycling does not lead to an agreement with the 
observed GSMF. A model that resembled vzw, including re- 
cycling, up to Ms ~ 10^1 Mq, but suppressed wind recycling 
with increasing effectiveness at higher masses would better 
match observations. 



Figure 4. The fractional stellar mass of central galaxies assem- 
bled via the different modes as a function of halo mass. Coloured 
lines are medians of the fractional stellar mass growth correspond- 
ing to each mode. Cold mode dominates the no-wind simulation 
(top panel) at all halo masses, always exceeding 75% of the to- 
tal mass growth. Cold mode exceeds hot mode in every feedback 
model (bottom three panels), but wind mode becomes the dom- 
inant growth mode at high Mj, in every case. We indicate the 
relation between halo and stellar mass by grey vertical lines indi- 
cating the median Mji for the given Ms listed next to each grey 
line. 



ends of the range plotted, similar to the findings of K09b. 
The vzw model traces this curve up to Mnw ~ 7 X fQi" Mq, 
just as it agrees with the observed GSMF up to this mass 
scale in Figure (5] However, none of our models follows the 
black curve over the full mass range, just as no model re- 
produces the observed GSMF at all masses. 

To compare the wind models to each other, we first con- 
sider the thin lines in Figure [6] which show /supp when wind 
recycling is explicitly removed - i.e. any particle that was 
ever in a wind is not counted when computing the galaxies' 
z = masses. The mass suppression in the cw and sw mod- 



4 DISCUSSION 

4.1 How Feedback Suppresses Star Formation 

As noted in our discussion of Figure[6l in the absence of wind 
recycling, stellar mass suppression in the sw and cw mod- 
els is considerably stronger than the naively expected factor 
of 1 -I- = 3. One may expect ejective feedback to reduce 
star formation by two thirds when — 2, because two Mq of 
mass is ejected for every one Mq of stars formed. Part of the 
difference reflects the impact of stellar mass recycling (dis- 
tinct from wind recycling). With our Chabrier IMF, 18% of 
the mass that forms stars is promptly recycled to the ISM as 
supernova and massive star ejecta, and over 10 Gyr this frac- 
tion rises to ~ 50%, mostly from AGB winds. In a no- wind 
model, this recycled gas generally forms stars, increasing the 
galaxy stellar mass. However, in a constant-wind model, a 
fraction of the initial gas supply is ejected and hence 
never contributes any recycled stellar mass, and of the 
gas that is returned to the ISM by evolved stars is ejected 
rather than forming stars. With stellar recycling fractions of 
50%, this effect alone can raise the mass suppression factor 
from 3 to 5. 
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Figure 5. The z = galactic stellar mass functions (GSMFs) for our four simulations including all accretion (thick l ines) and without 
wind mode (thin lines); the no wind (upper left panel) simulation, by definition, has no wind mode. Also shown is the I Bell et al.l ||2003| ) 
g-band derived GSMF (black) scaled to a Chabrier IMF for comparison. 



Furthermore, winds inject energy and momentum into 
their surroundings and thereby suppress accretion, provid- 
ing preventive feedback in addition to ejective feedback. This 
preventive feedback has the greatest impact on galaxies be- 
low M* . To investigate this, we ran a cw simulation where 
we did not allow winds to interact hydrodynamically at all. 
The stellar masses increased by a factor in excess of 2 for 
the smallest galax;ies we resolve. The energy added by these 
strong winds to the surrounding gas via hydrodynamic in- 
teractions is evidently playing a key role in suppressing gas 
accretion. This local "pre-heating" prevents IGM material 
from joining the filaments through which accretion is effi- 
cient. However, wind pre-heating is far weaker in the more 
moderate wind prescriptions. With only a quarter as much 
energy imparted to the surrounding gas, the sw winds result 
in galaxies twice as massive below Mnw = 10^^ Mq com- 
pared to the cw prescription (cf. the thin orange and cyan 
lines in Figure [6|. The vzw winds have even lower wind en- 
ergies at the low-mass end. Nevertheless, in all wind models 
cold mode SFRs decline faster at late times than expected 
from ejective suppression alone, likely owing to more dilute 
gas that can more easily be affected by energy input. 

The vzw case also behaves somewhat differently than 
simple expectations. Outflow suppression does not rise as 
significantly as one would expect toward higher masses (thin 
magenta line in Figure [6]), despite the r] oc scaling. If 
we assume u oc A^^''^ then M-"-^'^ (at a given redshift), 
which OD08 showed is a fair description for the vzw mass 



feedback rates. Over the range of Mnw = 10^°-^ - 10^^'^ M©, 
we would expect an increase of /supp by about a factor 4 — 5 
in the vzw suppression efficiency relative to the constant 
mass loading of the cw and sw models given how rj is calcu- 
lated. However the increase appears to be at most a factor of 
2, mainly because the more massive galaxies are the hierar- 
chical merger products of low-mass progenitors, which had 
higher rj before merging. This is the same argument K09a 
used to explain why massive galaxies are mainly assembled 
via cold mode - their progenitors are multiple small galaxies 
receiving gas via the cold mode channel. Thus, the higher rj 
of the progenitors balance the low rj of the merged galaxy 
at late times, reducing the differential mass loading effect. 

These results show that, even in the absence of the re- 
cycled wind mode, the suppression of stellar masses does not 
follow simple scalings with the mass loading factor Tj{A'Ib). 
Stellar mass recycling, hierarchical galaxy assembly, and, 
above all, the preventive feedback effects of winds (wind 
pre-heating), all act to complicate the mass dependence of 
feedback suppression. These results suggest caution when 
using analytic or semi-analytic models tied to the observed 
GSMF to infer the scaling of rj with galaxy mass. Further- 
more, we find that wind recycling changes the mass depen- 
dence of feedback suppression more strongly still, so model 
predictions will be sensitive to the assumptions (explicit or 
implicit) about the fate of gas after it has been ejected from 
the star-forming regions of its host galaxy. 
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Figure 6. The suppression factor /supp of galaxy stellar masses in 
feedback simulations relative to the stellar masses in the no-wind 
simulation. Galaxies are matched with the nw simulation by rank- 
ing them in order of decreasing stellar mass, then selecting galax- 
ies of the same rank. Thick curves show results for the cw (cyan), 
sw (orange), and vzw (magenta) models. Thin curves show the 
results when recycled wind mode is explicitly removed. The black 
solid curve shows the required /supp of no- wind stellar masses, 
Afnw , to match the masses of galaxies w i th the same comoving 
space density according from iBell et ah I ||2003| '). A model with 
■fsuDD following the black curve would reproduce the iBell et al] 
(120031 ) GSMF exactly. The dotted grey line shows the suppres- 
sion naively expected from ejective feedback in the cw and sw 
wind models when recycling is suppressed. 



4.2 Recycling and Faint-End Slopes of the GSMF 

Differential recycling and wind mode creates much more di- 
versity between the slopes of the mass functions (cf. thick 
coloured lines in Figure [5} than feedback suppression alone 
(thin lines). The sw model results in three distinct mass 
regimes, each with a different slope. At Ms < 10^" Mq 
the slope remains steep until the wind mode significantly 
softens the slope as it increases sharply between Ms ~ 
10^° - 10" Mq. Above Ms > lO" Mq the slope steep- 
ens again as efficient recycling causes it to approach the 
nw slope. This triple slope behaviour has also been iden- 
tified in the Overwhelmingly Large Simulations (OWLS) 
when using a constant wind-t ype model (Marcel Haas, pri- 
vate communication, see also ICrain et all (|2009l )). The cw 
mass function could be described in a similar way, although 
with higher v^i-nd the "flattening" from the wind mode is 
better described as a bump at high masses. Increasing Wwind 
by a factor of two dramatically alters the z = Q GSMF in 
two main ways: (i) preventive feedback suppresses by x2 
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Figure 7. The GSMFs of the vzw model from 2: = 3 in- 
cluding all accretion (thick lines) and without wind mode (thin 
lines). 



cold and hot mode stellar growth at most masses, and (ii) 
wind mode galaxy growth begins at ~ x8 higher mass. 

The vzw simu l ation manages to follow the faint-end 
slope of iBeU et al] (|2003| ) between Ms = 10^ - 10^° '' Mq 
fairly well. The flatter faint end is a direct result of the 
gradual differential recycled wind mode inherent in the vzw 
model. Above the threshold mass where recycling domi- 
nates, outflows behave as a series of fountains eventually 
becoming unable to escape their parent halo (i.e. the halo 
fountain phenomenon at > M*, OD08). Below the thresh- 
old mass where winds mainly escape into the IGM, wind 
mode contributes less to galaxies resulting in a very subtle 
steepening of the GSMF below lO''-'^ Mq, a milder version 
of the "triple-slope" behaviour in the sw model. Our reso- 
lution limit prevents us from concluding whether this simu- 
lated GSMF reproduces the "faint-e nd upturn" relativ e to a 
Schechter function, hinte d at by the Bell et al] (|2003l ) data 
and seen more clearly by iBaldrv et al.l l|2008l ): however we 
suggest the possibility this upturn could represent the tran- 
sition mass below which wind mode stops contributing and 
the slope steepens. 

Differential recycling also results in an evolution in the 
faint end of the GSMF, such that it is steeper at high red- 
shifts when recycling is a small effect compared to direct 
accretion. Figure [7] demonstrates this evolution for the vzw 
model. Feedback sets the slope of the GSMF at all epochs; 
the steeper slopes at high-z [z > 2) are set by the ejection 
of gas while the flatter slopes at late times are set by the 
recycling of gas. This evolution in the faint end has long 
been a challenge for hierarchical models, with only rela- 
tively ad hoc explanations proposed (see e.g. lKhochfar et"al] 
l2007t ). It is encouraging th at our vzw wind model can fit the 
observed high- ^: GSMFs in Oppenheimer. Dave, fc Finlatoij 
(|2009l . see also lDave. Finlator. fc Oppenheim er"2006''l . while 
concurrently showing a number of other successes over other 
wind models in reproducing key galaxy and IGM statistics 
(see §1). 



4.3 Is Wind Recycling a Numerical Artifact? 

The detailed hydrodynamic and radiative interactions that 
might drive a galaxy's ISM into outflows occur on scales 
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well below the resolution limit of any SPH simulations that 
model a representative cosmological volume. As discussed 
in !j2l Gadget-2 uses a "sub-grid" model to store thermal 
energy from stellar feedback in the hot phase of a two-phase 
ISM, which slows star formation, pressurises and stabilises 
the ISM, and improves numerical c onvergence by lower- 
ing th e sensitivity to mass resolution (jSpringel fc Hernguistl 
1200381 ). However, this ISM pressure does not, on its own, 
lead to galactic-scale o utflows, so Gadget -2, like some 
of its predecessors (e.g. iNavarro et al.l 1 19961 ). implements 
winds by imparting kinetic energy to individual particles 
(jSpringel fc Hernguistl l2003al ) . Winds are temporarily hy- 
drodynamically decoupled as described in !}2] to simulate 
the formation of hot chimneys extending out of a disk galaxy 
and to obtain resolution conve rgence in feedback simulations 
jSpringel fc Hernguistl boOSbl ). Without this decoupling, it 
is still possible to d rive winds that suppress star formation 
(ISchaye et all |2009| ) and enrich the IGM (|Wiersma et all 
l2009al ). but the results do not remain w e ll-con verged with 
resolution (e.g. iDalla Vecchia fc Schavd 120081 '). We have 
modified the Gadget-2 algorithm to implement the velocity 
scaling of momentum-driven winds (following OD08), but we 
have not altered the basic mechanism by which wind feed- 
back is implemented. Our exploration of feedback in this 
paper is limited to kinetic feedback only, and other met hods 
of feedback (e.g. SN blastwaves, IStinson et al.l I2OO6I I will 
result in different answers. 

By construction, our scheme achieves its original goal, 
ejecting gas from star-forming galaxies with mass-loading 
factors and wind velocities that are motivated by a combi- 
nation of theoretical and empirical arguments. Our results 
(and those of OD08) show that re-accretion of these ejected 
wind particles plays a major role in shaping the distribution 
of galaxy stellar masses and the correlations among halo 
mass, stellar mass, SFR, and redshift. 

Unfortunately, there are a variety of numerical effects 
that could cause our modelling of the post-ejection be- 
haviour of wind particles to be inaccurate. Most obviously, 
ejected wind particles initially have velocities and tempera- 
tures that are very different from those of their neighbours 
in the galaxy halo, while SPH is designed to represent fluid 
elements by groups of N ^ 32 particles. Second, SPH has 
known difficulties in treating two-phase media with sharp 
discontinuities in temperature and density, where gaps that 
develop around dense clumps tend to s uppress instabilitie s 
that might otherwise break them apart (|Agertz et al.ll2007l ). 
Third, the ram pressure on under-resolved clumps (or sin- 
gle particles) can be overestimated because their smooth- 
ing lengths must be large enough to enclose 32 neighbours, 
which may give them an artificially large cross-section. K09a 
suggest that this excess ram pressure causes some or all 
of the "cold drizzle" of T « 1 0^ K gas onto th e central 
galaxies of massive halos (see also lNaab et al.|[2007l '). Fourth, 
our calculation of meta l-line cooling is based on colhsional 
equilibrium tables while lWiersma et all (|2009al ') showed that 
photoionisation can significantly increase metal-line cooling 
times. In simulations without metal-line cooling, we find a 
significant decrease in wind recycling, and in overall gas ac- 
cretion for galaxies above M* , in gen eral agreement with the 
findings of IChoi fc Nagamind l|2009l ). Conversely, including 
the mixi ng of metals via inst abilities, conduction, or diffu- 
sion (e.g. IWadslev et al|[2008h would generally increase the 



overall accretion in hot halos, as me tal coolants spr e ad over 
more mass decrease cooling times. IWiersma et al.l l|2009bl ) 
demonstrated that using smoothed SPH metallicities to sim- 
ulate metal mixing can increase by 50% the number of 
baryons in stars compared to using single particle metal- 
licities. 

Wind recycling is unquestionably an important effect 
in our simulations, as we have shown in the preceding sec- 
tions. It is a physically plausible phenomenon that doubt- 
less occurs to some degree in the real Universe. However, 
given the numerical issues listed above, we cannot be sure 
that our simulations give an accurate quantitative account 
of the level of wind recycling that should arise in our 
adopted physical scenarios (cw, sw, and vzw). These nu- 
merical issues should be less severe for lower mass galax- 
ies without hot halos, therefore the predictions for re- 
cycling relating to the sub-M* galaxies and the GSMF 
faint-end slope should be more robust. It is worth re- 
iterating that simulations with the Gadget-2 wind al- 
gorithm and the vzw scalings have achieved a number 
of empirical successes, includi ng the z = 6 galaxy lu- 
minosity function (Dave . Finla tor. fc Oppenheimer 200^), 
the galaxy mass-metallicity relation ( Finlator fc Davd 
120081 ) ■ observations of IGM metal enrichment over a 
wide range of redshifts llOppenheimer fc Pavel l2006l . l2009l : 
lOppenheimer. Dave, fc Finlatoj|2009i r and the shape of the 
z = GSMF below M* (this paper). These agreements do 
not depend much on the recycling behaviour in high-mass 
galaxies. Conversely, wind recycling contributes to the exces- 
sive masses and SFRs of super-M* galaxies in high-mass ha- 
los, which is partially why the same simulations over-predict 
the global SFR density below z < 1 (Figure [1} and fail to 
reproduce the local ly observed colour bimodality of galaxies 
(|Gabor et al.ll2010l ). 

In comparing simulations with different mass resolu- 
tions (2 X 256^ particles in 16, 32, and 64/i~^Mpc boxes), 
OD08 found similar qualitative results for wind recycling but 
significant changes in tree at a fixed galaxy mass. Further in- 
vestigations of resolution and box size effects will be needed 
to characterise the numerical robustness of our quantita- 
tive results. In a different approach, we are using very high 
resolution simulations with idealised geometries to compare 
the behaviour of single-particle ejection to ejection of well 
resolved gas blobs (M. Peeples et al., in preparation). The 
strongest test of numerical robustness would be to carry out 
simulations that implement momentum-driven winds with a 
substantially different numerical algorithm; this is a longer 
term goal. At the same time, one can search for more tar- 
geted observational tests of our numerical predictions for 
galactic outfiows and wind recycling, a point we return to 
below. 



5 CONCLUSIONS 

We have examined the factors that govern the present- 
day galaxy stellar mass function in cosmological hydrody- 
namic simulations that incorporate several different galac- 
tic outfiow prescriptions. Our key result is the identifica- 
tion of a mode of accretion that we call recycled wind 
mode (or just wind mode), in which the accretion comes 
from material that was previously ejected from a galaxy. 
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This provides a third distinct accretion mode along with 
the " c old" and "hot" mo des d escrib ed by iKeres et alj (e.g. 
l2005h : lDekel fc Birnboiinl (e.g. I2OO6 I). At some level, the ex- 
istence of a recycled wind mode is obvious - material ejected 
from galaxies can be subsequently re-accreted - but the im- 
pact on our simulations is remarkably large. In our sim- 
ulation without winds, we recover the results of Keres et 
al. (2009a/b), in which cold mode accretion dominates the 
growth of galaxy masses. In simulations with winds, how- 
ever, wind recycling is the dominant mechanism of galaxy 
growth at 2: 1. 

Wind mode accretion occurs despite winds being ejected 
primarily at velocities exceeding the escape velocity of the 
galaxies' halos. Hydrodynamic processes are mostly respon- 
sible for slowing the winds (as also argued by OD08). Even 
when the outflow speed is proportional to galaxy velocity 
dispersion, winds recycle back into more massive galaxies 
faster because they live in environments of higher gas den- 
sity. The mass dependence of the wind mode we call differ- 
ential wind recycling. In all three of our wind models, the 
median recycling time declines at higher halo masses, but 
the rate of decline depends on the wind prescription. In the 
most massive halos, the recycling time is typically below 
1 Gyr, and > 90% of wind particles recycle without ever 
leaving their parent halos. 

The differential wind recycling relations ( i]3.2l) are mir- 
rored in the galaxy SFRs decomposed by accretion mode 
f ^3.3|) . Wind mode dominates the SFR of central galax;ies 
approximately when tree becomes less than the Hubble time, 
so that the outflow material is then accreted back into galax- 
ies at a rate exceeding the cold mode inflow. This transi- 
tion to wind mode dominance is reflected in the stellar mass 
growth of 2 = galaxies, as the galaxy stellar mass func- 
tions (GSMFs) from our three wind models show signiflcant 
differences that are directly traceable to the behaviour of 
the wind mode in each case. 

If we explicitly eliminate wind mode accretion in post- 
processing - by not counting recycled particles when com- 
puting galaxy stellar masses - then the z = GSMFs of our 
three wind models are surprisingly similar to each other, 
and completely unlike the observed GSMF. The suppres- 
sion relative to the no-wind simulation is much larger than 
the naively expected factor of (1 + rj), more so when high- 
velocity winds heat the surrounding gas and suppress sub- 
sequent accretion, providing preventive feedback in addition 
to the ejective feedback. When we include wind recycling (as 
our simulations do by default), then the GSMFs of the three 
wind models become quite different from one another, and 
the GSMF of the m omentum-driven feedback (vzw) model 
agrees best with the lBell et al.l (|2003l ) data for galaxy stel- 
lar masses 10^ Mq < Ms ^ 5 x 10^° Mq. The addition of 
wind mode flattens the faint-end mass slope from a = —1.92 
to -1.45 between lO^ '^ Mq Ms ^ 10^°'^ Mq. However, 
the predicted GSMFs greatly exceed the observed GSMF for 
Ms ^ 10^^ M0, approaching the prediction of the no- wind 
simulation, because material ejected in winds is quickly re- 
accreted in high mass halos. 

Reproducing the observed z — GSMF through wind 
mode alterations alone, at least in the vzw model, would 
require strong suppression of wind recycling at super-Af* 
masses (Mh ^ 2 x 10^^ M©), but minimal suppression of 
wind recycling at lower masses. It is conceivable that AGN 



feedback could provide a mechanism for such mass depen- 
dent suppression, just as it has been inv oked to suppress hot 
mode accretion in high mass halos (e.g.. ICroton et al]|2006l : 
ICattaneo et al.i .2009 ). However other mechanisms including 
virial shock heating (|Birnboim et al.ll2007l : lNaab et al.ll2007l ') 
and gravitat ional heating by infalling satellites and cold 
gas c lumps jPekel fc BirnboimI I2OO8I : iKhochfar fc OstrikeJ 
'2008) may suppress star formation without AGN, although 
the resolution in our cosmological-scale simulations may be 
insufflcien t to resolve su c h pro cesses. Using semi-empirical 
modeling iHopkins et alj l|2008l ) demonstrated that if star 
formation is efficiently quenched after major mergers of gas 
rich disks then the fraction of passive galaxies as a function 
of stel lar mass is well matched to observations. iGabor et al.l 
(|2010t ) showed quenching mechanisms applied post-run to 
vzw simulations can produce a red sequence of galaxies 
without dramatically changing galaxies below L* , although 
in detail none of the mechanisms tested there can exactly 
reproduce both observed luminosity functions and color- 
magnitude diagrams. Even if all star formation is quenched 
at M* and above, this cannot account for the factor of 
~ x5 greater SFR density than observed. Furthermore, 
iFirmani et al.l (2009) demonstrated recycling preferentially 
increases the specific SFR (SFR/Ma) in massive galaxies, 
which is in contradiction with observations showing the op- 
posite trend. Simulations including self-consistent mecha- 
nisms to quench massive galaxy star formation may yield 
unexpected consequences for the star-forming sequence of 
galaxies. 

We note that the impact of wind mode is particularly 
strong at late cosmic epochs, while at earlier times {z > 2) 
the GSMF faint-end slope is set more by ejection. In our 
favoured wind model this causes the faint end to become 
shallower with time, which is qualitatively in agreement with 
observations. We leave a more detailed examination of such 
evolutionary trends for future work. 

While wind recycling is an important effect in our sim- 
ulations, we discussed (in i]4.3p several reasons that our nu- 
merical treatment of winds and wind recycling could be in- 
accurate. Assessment of our predictions will require further 
testing of our sensitivity to numerical resolution, a study 
of the performance of the Gadget-2 wind implementation 
(particularly the single-particle ejection scheme), and con- 
sideration of alternative treatments of metal-line cooling. 
We can also look for direct observational tests of our nu- 
merical predictions of wind recycling and the distribution 
of ejected gas in galactic halos. One indirect line of evi- 
dence comes from th e galaxy mass-metallicity relation (e.g. 
iTremonti et al.ll2004h . which implies a steep dependence of 
mass-loading factor on halo virial velocity, roughly 77 oc v~j^ 
(M. Peeples & F. Shankar, in preparation; R. Wechsler, pri- 
vate communication). This dependence, steeper than that 
expected for momentum- or energy-conserving winds, could 
potentially be explained by a combination of the scaling 
for momentum-driven winds with the mass-dependent recy- 
cling rates that cause more massive galaxies to re-accrete 
much more of the enriched gas they once ejected. More 
direct tests of the simulation's predictions can come from 
quasar absorption studies that probe metals and ionization 
states in the environments around galaxies. The advent of 
the Cosmic Origins Spectrograph on Hubble Space Telescope 
allows much more sensitive searches for low-redshift metals. 
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including tests of the lOppenheimer fc Davd (|2009l ) predic- 
tion that most metals traced by O vi reside near galaxies 
and frequently recycle. 

Modellers have long argued that galactic outflows are a 
crucial element of galaxy formation, necessary to understand 
the shape of the galaxy luminosity function and the over- 
all low efficien cy with which cosmic baryons are converted 
into stars ( e.g. White & Frcnk 1991: Kauffman n et al1ll993l : 
ICole et all |l99 4; S omerville fc Pri mack 1999,)). Hydrody- 
namic simulations that explicitly incorporate winds into 
them are more successful than those without winds at re- 
producing other aspects of galaxy evolution and IGM en- 
richment. Those with moderate outflows from the majority 
of galaxies both effectively suppress star formation and en- 
rich the n earby IGM where most metals are observ ed to 
reside (e.g. IStocke et al.ll2006l : fWakker fc Savagdl2009D . Our 
results show that these outflows are by no means the end of 
the story - to understand how galaxies grow, one must un- 
derstand, theoretically and observationally, the fate of the 
gas that they expel. 
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